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Abstract. Curved broad-band spectral distributions of non-thermal sources like blazars are described well by a log-parabolic 
law where the second degree term measures the curvature. Log-parabolic energy spectra can be obtained for relativistic electrons 
by means of a statistical acceleration mechanism whose probability of acceleration depends on energy. In this paper we compute 
the spectra radiated by an electron population via synchrotron and Synchro-Self Compton processes to derive the relations 
between the log-parabolic parameters. These spectra were obtained by means of an accurate numerical code that takes the 
proper spectral distributions for single particle emission into account. We found that the ratio between the curvature parameters 
of the synchrotron spectrum to that of the electrons is equal to ~0.2 instead of 0.25, the value foreseen in the (5-approximation. 
Inverse Compton spectra are also intrinsically curved and can be approximated by a log-parabola only in limited ranges. The 
curvature parameter, estimated around the SED peak, may vary from a lower value than that of the synchrotron spectrum up to 
that of emitting electrons depending on whether the scattering is in the Thomson or in the Klein-Nishina regime. We applied 
this analysis to computing the synchro-self Compton emission from the BL Lac object Mkn 501 during the large flare of April 
1997. We fit simultaneous BeppoSAX and CAT data and reproduced intensities and spectral curvatures of both components 
with good accuracy. In particular, the large curvature observed in the TeV range was found to be mainly intrinsic, and therefore 
did not require a large pair production absorption against the extragalactic background. We regard this finding as an indication 
that the Universe is more transparent at these energies than previously assumed by several models found in the literature. This 
conclusion is supported by recent detection of two relatively high redshift blazars with H.E.S.S. 

Key words, radiation mechanisms: non-thermal - galaxies: active - galaxies: BL Lacertae objects. X-rays: galaxies: individual: 
Mkn 501 

1 . Introduction quency decades. A very simple and successful analytical func- 

tion that can model the shape of these broad peaks is a parabola 
The spectral energy distributions (SEDs) of blazars, and par- ^^e logarithms of the variables (hereafter log-parabola). This 
ticularly of BL Lac objects, are characterised by a double- f^^^^^^^ has only three spectral parameters, and was used by 
peak structure. The peak at lower energies is generally ex- Landau et al. (1986) to fit the broad band spectra of some 
plained by synchrotron radiation from relativistic electrons in a ^^-^^^^ Lac objects from radio to UV. Sambruna, Maraschi 
jet that IS closely aligned to the line of sight, while the high ^ ^rry (1996) also used a log-parabola to fit blazars' SED 
frequency bump is produced by inverse Compton scattering ^^^^^^^^ frequencies and luminosities without, how- 
from the same electron population interacting either with the ^^^^^ attributing a physical meaning to this shape. In two pre- 
synchrotron photons (SSC, Synchrotron-Self Compton mod- ^-^^^ ^^^^^^ (Massaro et al. 2004a, 2004b, hereafter Paper I 
els, see e.g. Marscher & Gear, 1985) or with other photons that p^p^^ respectively) we used the log-parabolic model to 

originated in the local environment (ERC, External Radiation j^e BeppoSAX broad band X-ray spectra of the BL Lac ob- 
Compton models, see e.g. Sikora, Begelman & Rees, 1994). -^^^^ 42I and Mkn 501 and studied the relations between 
For the so-called Low-energy peaked BL Lac (LBL) ob- the spectral parameters and the luminosity. We also showed 
jects, the frequency of the first peak is in the Infrared-Optical that, under simple approximations, a log-parabolic synchrotron 
region, while it is in the UV-X ray range for the High-energy spectrum can be obtained by a relativistic electron population 
peaked BL Lac (HBL) (Padovani and Giommi 1995). The having a similar energy distribution, and then derived the main 
shape of these bumps is characterised in the Log(vFy) vs. Log v relations between their parameters. Furthermore, in Paper I we 
plot by a rather smooth curvature extending through several fre- proposed a simple explanation for the log-parabolic energy dis- 

tribution of the electrons as resulting from a statistical mecha- 
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nism where the acceleration probability decreases with the par- 
ticle energy. 

Curved spectra in non-thermal sources have already been 
studied in the past and generally were related to the radia- 
tive ageing of the emitting electrons that have a single power 
law injection spectrum. The well-known solution by Kardashev 
(1962) predicts a change in the electron spectral index by a 
unity (0.5 for the synchrotron spectrum) around a break en- 
ergy, whose value decreases with time. Probably for this rea- 
son, curved spectra have been modelled by means of rather 
complex functions like a double broken power law (Kataoka 
et al. 1999) or a continuous combination of two power laws af- 
ter releasing the condition on the difference of spectral indices 
(see the discussion in Fossati et al. 2000). A similar spectrum 
is generally used in the spectral modelling of gamma-ray bursts 
(Band et al. 1993). On the same basis, Sohn et al. (2003) more 
recently introduced a spectral-curvature parameter, defined as 
the difference between two spectral indices in different fre- 
quency intervals, to study the evolution of synchrotron sources. 
As we will show later, this parameter is strongly dependent on 
the chosen frequency range and cannot be univocally related to 
the actual curved shapes. 

In this paper we stress the advantages of using a log- 
parabolic law to model curved spectra over broad frequency 
ranges and present the results of an accurate study of the spec- 
tral distributions of the synchrotron radiation (hereafter SR) 
and of the inverse Compton (hereafter IC) scattered photons 
by a population of relativistic electrons with such energy spec- 
trum. We study, in particular, the relations between the parame- 
ters of SR spectra and those of the electrons and show that they 
are generally approximated well by simple power behaviours 
and that they provide useful information on the electron spec- 
trum. The relations with the IC spectrum are not so simple 
depending upon whether the majority of scatterings is in the 
Thomson or in the Klein-Nishina regime. In any case, our cal- 
culations indicate a good way to properly model the observed 
spectra. We also consider the case in which the probability of 
statistical acceleration is constant below a critical energy, thus 
producing a power law spectrum that turns into a log-parabola 
above this energy. 

In the last part of the paper, we apply our model to the large 
flare from Mkn 501 observed in April 1997. This is a unique 
event in which the SR and IC components were observed si- 
multaneously over energy bands that were broad enough to 
study the spectral curvature. We use the BeppoSAX data al- 
ready analysed in Paper II together with those in the TeV band 
obtained by the CAT experiment (Djannati-Atai et al. 1999). 
Our approach of studying both the SR and IC spectral curva- 
tures is useful for evaluating the relevance of the intergalactic 
absorption due to pair production from photon-photon interac- 
tions. In particular, we discuss the possibility that the observed 
TeV curvature is mainly intrinsic, as it implies rather low ab- 
sorption and the possibility of detecting sources at larger dis- 
tances than usually estimated. 



2. Properties of the log-parabolic spectral 
distribution 

The log-parabolic model is one of the simplest ways to rep- 
resent curved spectra when these show mild and nearly sym- 
metric curvature around the maximum instead of a sharp high- 
energy cut-off like that of an exponential. In the following, we 
sunmiarise the main properties of log-parabolic spectra, while 
for details we refer to Papers I and II. This law has only one 
more parameter than a simple power law and can be written as: 

F(E) = KiE/Ei)-^" ^ ^ (1) 

The energy dependent photon index is 

a(E) = a + 2b LogiE/Ei) . (2) 

With Vp - Ep/h we also indicate the peak frequency of the 
SED where Ep is related to the spectral parameters of Eq. (1) 
as follows: 

Ep = El lO*^-"'/^* (3) 
and 

VpFiVp) = KEiEpiEp/Er)-"'^ = KE\ K^-^fl^b . (4) 

It is useful to note that the above formulae must be used 
when working with photon spectra, as is usual in X and y-ray 
astronomy. Energy spectra are obtained by multiplying Eq.(l) 
by E, and the corresponding spectral index is therefore in- 
creased by unity. In the following, we indicate the energy in- 
dex with a<, = a - 1 to avoid confusion. In this case the peak 
frequency becomes: 

Vp = VI io(i-««)/2* (3'). 

An advantage of the log-parabolic model compared to other 
curved spectral laws is that the curvature around the peak is 
characterised only by the parameter b, while in other mod- 
els, like the continuous combination of two power laws, it is 
a function of more parameters. Note also that the estimate 
of b is independent of the particular spectral representation 
adopted for the data. The same property does not hold for 
the spectral-curvature parameter (SCP) defined by Sohn et al. 
(2003). Unlike b, this quantity depends upon the energy (fre- 
quency) interval chosen to estimate the spectral indices at high 
Oh and low ai frequencies. The relation between b and SCP can 
be obtained by evaluating a in two energy ranges E^ (high) and 
El (low) from Eq.(2): 

SCP = = ^ LogjEh/Ei) 

ah + ai a + b Log{EhEilE\) 

For flat spectrum sources with a close to zero, SCP becomes 
practically independent of b. 

One limit of the log-parabolic model, however, is that it can 
represent only distributions symmetrically decreasing with re- 
spect to the peak frequency. However, it is not difficult to mod- 
ify Eq.(l) to take into account a possible asymmetry with re- 
spect to Ep. For instance, one can use two different values of b 
for energies lower and higher than Ep, where the two branches 
continuously match. An interesting possibihty, which can be 
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simply explained in terms of the energy dependence of the par- 
ticle acceleration probability (see next section), is that the low 
energy segment of the spectrum follows a single power law 
with photon index ao and that the log-parabolic bending be- 
comes apparent only above a critical value Ec- This behaviour 
can be described by the following model that takes into account 
the continuity conditions on the flux and on a(E) at Ec: 



(6) 



FiE) = KiE/Eiy"" , E<Ec 
F{E) = KiE/Ecf^^'-'^'"^'^ X 

3. Statistical particle acceleration and 
log-parabolic spectra 

In Paper I we showed that a log-parabolic distribution is not 
only a simple mathematical tool for spectral modelling, but that 
it also relates to the physics of the statistical acceleration pro- 
cess under some simple hypotheses. In the following, we ex- 
tend the considerations of Paper I. The results, useful to derive 
some relations between the acceleration parameters, should not 
be considered as a complete theory of energy-dependent accel- 
eration. To confirm the main findings of this approach, more 
detailed analytical and numerical calculations are necessary, 
but such a detailed theoretical approach is beyond the purpose 
of the present paper. 

3.1. Energy distribution of accelerated particles 

The energy spectrum of particles accelerated by some statisti- 
cal mechanism, e.g. a shock wave or a strong perturbation mov- 
ing down a jet, is given by a power law (Bell 1978, Blandford 
& Ostriker 1978, Michel 1981). In Paper I we showed that a 
log-parabolic energy spectrum is obtained when the condition 
that the acceleration probability p is independent of energy is 
released and that its value at the step / satisfies a power law 
relation as: 



Pi^g/rl , (i = o,i,2,...) 



(7) 



where g and q are positive constants. Such a situation can occur, 
for instance, when particles are confined by a magnetic field 
with a confinement efficiency that decreases for an increasing 
gyration radius. In Paper I (Sect. 6. 1) we found that the integral 
energy distribution of accelerated particles is given by a log- 
parabolic law: 



NO y) = A?o(r/ro)"''"'^''^''*'''^''°" 



with 



Log E 2 

q 



2 Log e 



(8) 

(9) 
(10) 



where yo is the minimum Lorentz factor in Eq.(7). As already 
discussed in Paper I the differential energy distribution A^(7') is 
not an exact log-parabolic law, but differs from it only by a log- 
arithmic term and can be approximated well by a log-parabola 




14 16 
Log(v) (Hz) 



Fig. 1. Synchrotron spectra emitted by an electron population 
with log-parabolic (dashed lines) and power-law log-parabolic 
energy distributions (Eq. 11) (solid lines), computed for r val- 
ues in the interval 0.50-1.20. Spectra were shifted on the verti- 
cal scale to avoid confusion. 

with a very similar curvature over an energy range of several 
decades. 

Shock wave acceleration is not the only statistical mecha- 
nism active in non-thermal sources. For instance, electron ac- 
celeration can occur in magnetohydrodynamical turbulence in 
which regions of magnetic field reconnection can develop in 
a very stochastic way. Recently, Nodes et al. (2004) presented 
the results of the numerical simulations of a relativistic particle 
acceleration in a three-dimensional turbulent electromagnetic 
field configuration, also taking their SR into account. These 
authors found energy spectral distributions that where signif- 
icantly flatter than s - 2 and, in a few cases, characterised by a 
steepening spectral index at high energies. We verified that over 
sufficiently wide energy ranges, the spectra given by Nodes et 
al. (2004) are represented well by a log-parabolic law or by 
a combination of a power law and a parabola. The resulting 
curvature parameters are generally small, but probably the cur- 
vature depends on the distribution and size of the acceleration 
regions, so one can expect that, under different assumptions, it 
could be higher Nodes et al. (2004) computed the emerging 
SR spectra, which show an appreciable curvature over a few 
decade frequency range. 

Statistical acceleration is not the only way to obtain curved 
electron spectra. Energy distributions showing a rather mild 
curvature have also been obtained in blazar physics. To model 
the SED of MeV blazars, Sikora et al. (2002) assumed that 
electrons are accelerated via a two-step process with a broken 
power-law energy distribution as injection. When the cooling 
effects are taken into account, the resulting electron spectrum 
(see, for instance. Fig. 6 in their paper) can be described well 
by a log-parabola over a range that is wider than three decades, 
as we verified (in this case we found r=0.44). 
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The assumption of Eq. (7) about the energy dependence of 
the acceleration probabiHty can be modified to take other es- 
cape processes into account. For instance, one can assume that 
the acceleration probability is constant for low energies and 
that it begins to decrease above a critical Lorentz factor ja (see, 
for instance Eq. (23) in Paper I). An approximate expression for 
the energy distribution of accelerated particles under this con- 
dition for y <s^ ja should follow a power law with spectral index 
So - -{Log g)/{Log s), while for y » ja it will approximate a 
log-parabolic spectrum like Eq.(8): 



N{y) = No (r/ro) r < ro 

N{y) = No (r/ro)"*-'"^'' '-"^^^'T"^'' y > 



ro 



(11) 



where for the sake of simplicity we take yo-ya- In addition to 
the power law and log-parabolic segments, we can also expect 
a sharp high-energy cutoff that appears when losses make the 
acceleration process highly inefficient. 
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4. Spectral properties of the emitted radiation 

4.1. Synchrotron spectrum 

It is important to know the relations between the spectral pa- 
rameters a and b of the emitted radiation and those of the elec- 
tron population, namely s and r. The spectral distribution of 
the SR by relativistic electrons having a log-parabolic energy 
distribution cannot be computed analitically. For our purpose, 
however, the relations between the spectral parameters can be 
derived under the usual (5-approximation and the assumption 
that the electrons are isotropically distributed in a homoge- 
neous randomly oriented magnetic field with an average inten- 
sity B 



js(v) 



P(v(y))N(y) dy 



(12) 



where js is the synchrotron emissivity, N{y) the electron den- 
sity, and the power radiated by a single particle is 



P(v) - — 
9m?-c 

and 

vs = 0.29 



y^B^ 6{v-vs) 



3e_ 

Annie 



y^B 



(13) 



(14) 



is the SR peak frequency for a single particle. We then obtain 
(see Paper I): 



Ps(v)^NoB\v/vo) 
with 

a = (s-l)/2 
b = r/4 . 



-{a+bLog(v/vo)) 



(15) 



(16) 



Note that the parameter a, as defined above, does not coin- 
cide with the one defined in Eq. (2), because the former is the 
energy index at the frequency vo, whereas the latter is the pho- 
ton index at the energy E^, not necessarily corresponding to vq. 
It is not difficult to verify that these two parameters differ for an 



Fig. 2. The relation between the peak frequencies of the syn- 
chrotron SED and those of the energy spectrum of the emit- 
ting electrons computed using the SFD (Eq. 17) for the LP and 
LPPL cases (see text). The abscissa is the inverse of the elec- 
tron spectrum curvature to show the linear trend expected from 
Eq.(18). The dashed line is the best fit computed for 1/r >0.33. 
Note the small deviations of LPPL points from the best fit line 
at low values of 1 /r. 



additional constant, whose value depends on vq and vi - Ei/h. 
In Paper I we noted that one can expect a and b to show a linear 
correlation, which was indeed found for Mkn 421. We stress 
that such correlation is not a general result that can be applied 
to any sample of sources. It also depends upon the values of 
q, g, and yo, which can be different from source to source, and 
the dispersions of their distribution functions can introduce a 
significant amount of scatter in the data that could destroy the 
correlation. Reasonably, one has to expect that the linear cor- 
relation between a and b can be found only when the physi- 
cal conditions in the acceleration region are quite similar. The 
absence of correlation found by Perlman et al. (2005) in their 
sample, therefore, does not imply that curved spectra in blazars 
are not originated by energy-dependent statistical acceleration. 

The relation between r and b given above (Eq.l6) is not 
exact; one can expect that a precise calculation of the spectral 
curvature of the emitted radiation must give a b value smaller 
than in the 6 approximation and that, for high values of r, b can 
be greater than unity (depending on the frequency interval used 
around the peak in which it is estimated) because of the expo- 
nential cutoff by the SR spectrum radiated by a single particle. 

To compute the relations between the spectral parameters 
numerically and to study the spectral evolution of the emit- 
ting particles, as reported in the next subsection, we applied a 
time-dependent numerical code (Tramacere & Tosti, 2003; see 
Tramacere, 2002, for a detailed description of the code). This 
code uses a numerical integration routine that applies a modi- 
fied Simpson rule ensuring good convergence in short comput- 
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Fig. 3. The relation between the log-parabolic curvature pa- 
rameters of the electron energy distribution r and those of the 
synchrotron radiation computed applying the 6 approximation 
(crosses) and the SFD. The two LP point sets correspond to log- 
parabolic electron spectra with b evaluated by a best fit around 
the SED peak (filled circle) and on the high frequency branch 
(diamonds). LPPL points (open squares) correspond to a com- 
bined log-parabolic power law distribution. Dashed lines are 
linear best fits used for the calculation of the /j/r ratio. 

ing times. In particular, this code includes an accurate evalua- 
tion of the synchrotron frequency distribution (hereafter SFD): 



P(y) = 




where K^i^^x) is a modified Bessel function of the second 
type and Vc the synchrotron critical frequency (Blumenthal and 
Gould 1970, Rybicki & Lightman 1979). Note that while the 
maximum of P{v) is at 0.29 Vc, the peak in the corresponding 
SED vP{v) is at the frequency Vps =1.33 v^. 

We calculated the SR spectra and SEDs from Eq.(12) using 
a log-parabolic distribution without (LP) and with a low-energy 
power law branch (LPPL) (Eq. 11) and the frequency distribu- 
tion of Eq.(17). All the calculations in this and in the following 
subsections are performed in the co-moving frame, so relativis- 
tic beaming effects were not included. 

Figure 1 shows some examples of the resulting SEDs. 
In these computations we assumed a uniform magnetic field 
fi=0.14 G, and the spectral parameters for the electrons were 
70=10^, s=so=l-2 (Eq. 1 1), while r values were increased from 
0.5 to 1.20 by steps of 0.05. The electron density and the size 
of the emission region were taken to be equal to the typical 
values used in blazar modelling. In these conditions the syn- 
chrotron self-absorption is important at frequencies lower than 
~10'° Hz. Note in Fig. 1 that at low frequencies both spectral 
sets tend to a power law: in the LP case the spectral index cor- 
responds to the SFD asymptotic behaviour (v'^^), while in the 
other case it has the energy spectral index a^o - {sq - l)/2. 



At high frequencies the log parabolic behaviour is very evident 
with a peak frequency decreasing for higher values of r, as ex- 
pected from Eq.(3) since b is also increasing with r. 

We first verified that the peak frequency in the SED of 
the SR component Vps is proportional to the Lorentz factor jp 
at which the electron distribution A^(7') has its maximum. 
Applying the same formulae of Sect. 2, we have 

from which we have: 

^ = ^ 10'/^ (18) 

n To 

where we assumed vq = vi for the sake of simplicity. In Fig. 2 
we plotted the resulting Vps /j^ vs 1/r for both LP and LPPL 
cases: points of the two sets are coincident in the whole range 
of interest and only for 1/r < 0.3 a small deviation from the 
linear trend of the LPPL model is evident. The slope of the best 
fit line (for 1/r > 0.3) is equal to 0.97, very close to unity as in 
Eq.(18). Note also that the linear extrapolation to large values 
of r (1/r ^ 0) gives for the Vps /jp ratio a value very close to 
1.33 Vc, consistent with the theory expectation. 

We also used the spectra of Fig. 1 to study the relation be- 
tween r and b under different conditions. The estimate of b is 
not simple because the radiation spectra do not have an ex- 
act log-parabolic shape. One can expect that this estimate is 
more accurate considering only the SED branches at frequen- 
cies higher than the peak: a quadratic best fit of the spectra of 
Fig. 1 confirms the linear relation between b and r, but the coef- 
ficient is equal to 0.22, slightly lower than the value found un- 
der the 6 approximation (Eq. 16). We also computed b consider- 
ing a frequency interval centred at the peak frequency and span- 
ning approximately two decades. The results are shown in Fig. 
3. The relation between b and r was linear with very good accu- 
racy, but the coefficient decreased to 0.18, as expected because 
of the smaller curvature of the low frequency branches. The 
last result does not change if the low energy electron spectrum 
follows a power law instead of a log-parabola (open squares in 
Fig. 3): the b/r ratio in this case was found equal to 0.19. We 
also verified our results by computing the spectra using the 6 
approximation and found /7/r=0.2508, in very good agreement 
with the theory. An approximate rule derived from this analysis 
indicates that it is convenient to take the value of the ratio equal 
to 0.2, within an accuracy of 10%. The resulting r estimate can 
then be used as the input value for a more refined modelling of 
the SED. 

4.2. Spectral evolution under synchrotron cooling 

We also studied the evolution of spectral curvature due only to 
synchrotron cooling. Of course, this study cannot be performed 
in general because the solution of the time-dependent continu- 
ity equation is a function of the assumed time evolution of the 
injection rate of emitting electrons and of several other param- 
eters like the mean residence time of particles inside the emit- 
ting volume. We studied, therefore, the simplest case of a pure 
radiative SR cooling in a homogeneous and steady magnetic 
field without leakage and with a log-parabolic initial spectrum. 
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Fig. 4. Time evolution of the SR spectrum (upper panel) and 
of the corresponding spectrum of electrons, multiplied by the 
square of the Lorentz, under SR losses. The injection spectrum 
was a log-parabola with s = 1.2 and r - 0.7, corresponding to 
a peak energy of ~ 10^. Spectra are plotted for equal time steps 
of ~0.4 cooling times. Note the decrease in the peak energy and 
the deviation from a log-parabola due to difference in cooling 
time with y. 




Fig. 5. Time evolution of the log-parabolic curvature parame- 
ters r of the electron spectrum (upper line) and b of the SR SED 
(lower line) computed by means of a best fit around the peak. 
Solid lines are linear best fits to the computed values. 



The time evolution of SR and electron spectra are shown in 
the two panels of Fig. 4, where the decrease in both the peak 
frequency and intensity with time is very evident. Although ra- 
diation losses modify the spectral shape, it remains approxi- 
mately log-parabolic, though with larger curvature parameters. 
These variations are shown in Fig. 5, where we plotted the r 
and b values, measured in energy/frequency intervals centered 
approximately at the distribution peaks, as a function of time, 
measured in units of the cooling time for the peak energy of the 
initial electron spectrum, taken as y=10^ and corresponding to 
~ 5x10^ s with a magnetic field B = 0.1 G. It is worth noting 
that the change in time of r and b is described well by linear re- 
lations. The observation of a decrease in Vps associated with an 
increasing curvature during the dimming phase of an outburst 
is an indication of a radiative cooling evolution. Our calcula- 
tions also show that the changes of SR peak frequency and b 
are correlated according to a power law as shown in the plot of 
Fig. 6. 

We stress, however, that the estimate of b may not be sim- 
ple when working with observational data because the position 
of the peak may be far from the centre of the detector frequency 
range, and the uncertainty on the interstellar and intrinsic ab- 
sorption may be non-negligible, particularly in the soft X-ray 
band. 

4.3. Synchrotron Self-Compton spectrum 

We computed the spectral distribution of the SSC emission 
from a relativistic electron population with a log-parabolic en- 
ergy spectrum in the framework of a homogeneous one-zone 
model in which the seed photons are isotropic in the bulk frame 
of the electrons. Our goal is to estimate the relations between 
the main spectral parameters: the peak frequencies of the SR 
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Fig. 7. The SR and SSC spectral energy distributions emitted 
by an electron population with power-law log-parabolic distri- 
butions (Eq. 11), computed for r values in the interval 0.50- 
1.20. Upper panel: spectra for yo= 10^; lower panel: spectra for 
ro=2xlOl 

and IC components and their curvatures indicated in the fol- 
lowing by bs and be, respectively. 

The IC emissivity is given by (Jones 1968, Blumenthal & 
Gould 1970, Band & Grindlay 1985) 



v')c J dvn{v) J dy(do-/dv')N(y) (19) 



jciv') = (h 



where n(v) is the density of SR photons, v' the frequency of 
scattered photons, and the differential Klein-Nishina cross sec- 
tion is: 
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3crj 



16/1- 



.i(l 



2t] In77-H(l -h 277)(1 -77)- 
[4(hv/mec'^)'yi]f 



1) 



1 -t- 4(hv/mi.c^)yrj 



V 



4y^v[l - (hv'/m^c^)] 



< 1 



4y2 



1<^<- 

V 1-1- Ay{hv I nieC^) 



(20) 



(21) 



(22) 



Figure 7 shows two sets of resulting spectral distributions 
of both SR and IC components for values of r ranging from 
0.5 to 1.2 and 70 equal to 10^ (upper panel) and 2x10^* (lower 
panel). The upper cut-off of the energy electron spectrum was 



Fig. 8. Dependence of the peak frequencies upon the particle 
characteristic energies in the SED of a single zone SSC model. 
Upper panel: frequency of SSC peak vs electron Lorentz fac- 
tors of the peaks for r values in the interval 0.50-1 .90, the three 
curves correspond to 70= lO'' (filled circles), 5x10^ (open cicr- 
cles), 10'* (crosses). Lower panel: ratio between the peak fre- 
quencies of IC and SR components for the same values of the 
spectral parameters. 

fixed at ye,=5xlO^. In the upper panel of Fig. 8, we plotted 
the peak frequency of the IC component Vpc vs yp for yo=10^, 
5x10^, lO'* and r ranging from 0.5 to 1.9. These curves are 
useful for understanding at which electron peak energy the IC 
emission is dominated by the Thomson or the Klein-Nishina 
regime. We stress that this does not imply that the largest con- 
tribution to the emission around Vpc is due to electrons exactly 
at yp. We computed these energies and found that, for the high- 
est value of 70, the difference is on the order of a few percent, 
while for the lowest 70 and r=1.7 the greatest contribution to 
the IC peak is produced by electrons having 7 ^ 6.3x10^, while 
yp is at 1.4x10"*. It is of course impossible to derive a sim- 
ple rule because of the number of parameters involved. Figure 
8 (upper panel) shows that as yp increases, the corresponding 
Vpc grows slower and slower indicating that the Klein-Nishina 
suppression becomes more efficient. A second way to show the 
relation between SR and IC components is given in the lower 
panel of Fig. 8, where we plotted the ratio Vpc/vps vs 7^: this 
ratio increases until the Thomson scattering is dominant and 
reaches the maximum at the transition to the Klein-Nishina 
regime. 

The IC spectra show an evident curvature that is not uni- 
vocally related to the r as in the SR case. Spectral curvature, 
in fact, depends on several parameters in a complex way. The 
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spectral shape is generally different from a log-parabola, and 
it can be approximated by this function only in a limited fre- 
quency interval. Consequently, the estimate of the curvature 
parameter depends on the postion and amplitude of this inter- 
val. In particular, the curvature depends on the intrinsic electron 
spectral curvature, the energy of SR photons, and on the energy 
of the electrons that mostly contribute to the emission in the 
selected interval. These energies, in fact, determine if the scat- 
terings happen in the Thomson limit or in the Klein-Nishina 
regime. We obtained spectra with a curvature less pronounced 
than SR only for 70=10^ (i.e. for dominant Thomson scatter- 
ing) and b evaluated in an interval around the peak. In all the 
other cases, higher b values of IC spectrum resulted. Figure 
9 shows some examples: the curvature parameter b of the IC 
spectrum, computed in a single zone SSC model, was evaluated 
in three adjacent frequency intervals, each having an amplitude 
of about a decade starting from the peak frequency. When the 
main contribution to IC emission comes from Thomson scat- 
terings, the curvature is closer to that of SR (solid line), and it 
approaches that of the electron spectrum (dashed-dotted line) 
when the fraction of interactions in the Klein-Nishina regime 
increases. The two panels correspond to yo values differing by 
a factor of 20. In the upper panel, IC emission is dominated by 
Thomson scattering in the two first intervals, while in the lower 
panel the curvature approaches that of the electrons because the 
majority of interactions are in the Klein-Nishina regime. This 
property is really useful because a simultaneous measure of the 
curvature parameters of SR and IC emissions can help to dis- 
criminate between the two regimes of Compton scattering and 
to constrain the spectral parameters of emitting electrons. 



5. The X-ray and TeV emission from l\fll<n 501 

In Paper II we presented the results of analysing all BeppoS AX 
observations of the nearby HBL object Mkn 501 using log- 
parabolic spectral models. In the case of this source, we did 
not detect significant variations of the X-ray flux during the 
individual pointings, but large luminosity changes were found 
between the three observations of April 1997 (on the 7, 1 1, and 
16) for which simultaneous TeV data, obtained with the CAT 
Cerenkov telescope, are available (Djannati-Atai et al. 1999). 
The TeV data showed a remarkable spectral curvature mod- 
elled by these authors by means of a log-parabolic law: best 
fit b values were higher than 0.4, while those measured in the 
X rays are in the range 0.12—0.17 (see Paper II). Krennich et 
al. (1999) observed Mkn 501 from February to June 1997 with 
the Whipple telescope and found a high curvature of the TeV 
spectra as well. 

In Sect. 4 we discussed the relations between the parame- 
ters of log-parabolic spectra for a SSC emission and showed 
that SSC emission model can give high frequency spectra with 
a b value ranging from that of the SR component to that of the 
electron spectrum depending on whether the scattering occurs 
in the Thomson or in the Klein-Nishina regime. From the cal- 
culations of Sect. 3, we know that an SR b value of about 0.15 
corresponds to r ^0.75, and therefore a curvature of the IC 
component as that measured in the TeV range can be expected. 




" I , I , I , I I L_ 

0.4 0.6 0.8 1 1.2 
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Fig. 9. The curvature parameter b of the IC spectrum for a sin- 
gle zone SSC model plotted against the electron spectral curva- 
ture r. Dashed-dotted lines are the curvature parameters for the 
electrons and solid lines for SR component. Dashed lines cor- 
respond to IC curvatures evaluated in three adjacent frequency 
intervals with an amplitude of about a decade selected start- 
ing from Vpc'- lowest lines correspond to the intervals closest to 
the peaks and the highest ones to those at the highest frequen- 
cies. The upper panel shows b values computed for yo^lO^ and 
lower panel those for yo=2xlO'*. Note that when the largest 
contribution to IC emission is from Thomson scatterings (up- 
per panel), curvatures are closer to those of SR, while in the 
Klein-Nishina regime (lower panel) b values result very close 
to r. 



In Paper II we also presented a two component-model to 
describe the spectral evolution of the April 1997 X-ray outburst 
of Mkn 501. We assumed that the components had the same 
curvature (b =0.18) and differ for the peak energy and their 
relative intensity. In this way we explained the nearly constant 
flux observed at lower energies than ~1 keV, while that around 
100 keV changed by more than one order of magnitude, and 
the slight decrease of b in the higher states was interpreted as 
the effect of the sum of two log-parabolic distributions peaked 
at different energies. 

An important effect to be considered when computing the 
spectral curvature in the TeV range is the absorption of the most 
energetic photons due to pair production against the infrared 
extragalactic background light (EEL). Dwek and Krennich 
(2005) have recently studied twelve EBL models and computed 
the corresponding intrinsic spectra of Mkn 501. These authors 
found that the majority of EBL models give intrinsic emitted 
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V (Hz) 

Fig. 10. Two spectral energy distributions of Mkn 501 during the low and high states observed on 7 and 16 April 1997, respec- 
tively. X-ray points are from Paper II, TeV points are simultaneous CAT data (Djannati-Atai et al. 1999) and soolid lines are 
the spectra computed in a 1-zone SSC model for the SR and IC components. In the upper panel IC, spectra have been absorbed 
(dashed lines) by interaction with infrared EBL photons according to the LLL model by Dwek and Krennich (2005). In the lower 
panel EBL absorption was neglected. 



spectra showing an unphysical rise at high TeV energies and 
that only the lowest intensity EBL models are acceptable. 

We applied our radiation code to searching for the geomet- 
rical and physical parameters of the emitting region that can re- 
produce the observed intensities and spectral shapes, also tak- 
ing EBL absorption into account. To compare the model with 
the data, we need to introduce a relativistic beaming factor <5, 
which does not affect the spectral shape, but which changes 
observed frequencies and fluxes by S/{1 + z) and [S/{1 + z)]^, 
respectively, where z - 0.034 is the source redshift. In the fol- 
lowing we assume 5 = 15, in agreement with other models for 
Mkn 501 (see e.g. Kataoka et al. 1999). 

We first considered a single-zone SSC model and tried to 
reproduce the two extreme states of the large outburst of April 
1997, i.e. the low state on April 7 and the high state on April 
16. 

The spectrum of the emitting electrons was chosen with 
a curvature parameter r-b/Q.2^ 0.75 following the relation 
found in Sect. 4.1, and the parameter s was fixed to 1.2. Note 
that this value of r is consistent with the general behaviour 
presented in Fig. 9. The TeV curvature measured around the 
peak was about 0.4 for the low state and 0.45 for the high state 



(Djannati-Atai et al. 1999) and, for yo on the order of 10"* and 
even higher, we expect a value of r in the range 0.5-0.8. The 
values of the other parameters, given in Table 1, were chosen 
to achieve a satisfactory agreement between the data and the 
model. It is interesting to note the behaviour of the two cor- 
related parameters v^^ (the peak frequency in the synchrotron 
SED) and jj, (the value of y of the electrons that contribute 
mostly to synchrotron peak emission). In fact the Vps grows 
during the flare by a factor of 1 1 compared to the initial value 
while, according to Eq. (14), jp changes as the square root of 
Vps growing rate. Spectral distributions of Fig. 10 show that 
the SR follows the X-ray data accurately in both the spectral 
shape and peak evolution. The TeV spectra in the two states 
were computed using the SSC code described in Sect. 4.2: in 
the upper panel of Fig. 10 we considered the EBL absorption 
for the LLL model by Dwek and Krennich (2005), while in 
the lower panel this absorption was neglected. In both cases 
TeV data were well fitted in both the peak position and spec- 
tral curvature. This suggests that this curvature could be partly 
intrinsic rather than produced entirely by EBL absorption. The 
LLL model could then be considered an upper limit to the ex- 
tragalactic background. 



10 



E. Massaro et al.: Log-parabolic spectra in blazars 111: SSC emission in the TeV band from Mkn 501 




Fig. 11. Two spectral energy distributions of Mkn 501 during the high states observed on 7, 11 April 1997 (upper panel) and 
7, 16 April 1997 (lower panel). X-ray points are from Paper II, and TeV points are simultaneous CAT data (Djannati-Atai et al. 
1999). Thin solid lines are the spectra computed in a 2-zone SSC model for the SR and IC components, dashed lines are the 
spectra of the high-energy flaring component, and the thick solid line is that of a slowly evolving component. 



Table 1. Model parameters for the 1997 April flare of Mkn 501 . 



Date 


s 


r 


N 


70 


7p 


V 




vLyiSR) 


VLy{IC) 


B 


6 








cm 


10* 


10= 


cm' 


lO-*" erg 


lO^'erg/s 


lO^'erg/s 


G 




SSC - 1 zone 
























Apr. 07 


1.12 


0.81 


4.25 


4.80 


1.67 


1.57X10-*' 


91.0 


8.70 


4.47 


0.10 


15 


Apr. 16 


1.18 


0.70 


30.0 


8.70 


3.36 


1.72x10-*' 


90.6 


34.60 


21.4 


0.14 


15 


SSC - 2 zones 
























Apr. 07 (Zl) 


1.14 


0.81 


4.6 


4.0 


1.35 


1.2x10-*' 


57.9 


8.71 


3.42 


0.14 


15 


Apr. 1 1 (Z2) 


1.2 


1.00 


1050 


14 


3.52 


9.2x10-*° 


0.156 


6.92 


1.44 


1.5 


15 


Apr. 16 (Z2) 


1.2 


1.33 


250 


20 


3.95 


9.2x10*'' 


6.12 


29.5 


14.1 


0.2 


15 



V is the volume of the emitting region; is the total energy of radiating electrons; vLy(SR,IC) are the SED peak values of SR and IC 
components, respectively. 



As proposed in Paper II we considered a two-zone SSC 
model. A first zone (Zl) is responsible for the "slowly" vari- 
able emission, whereas the second zone (Z2) is the source of 
the high-energy flare: the observed fluxes are then the sum of 
these two contributions. We used the April 7 emission as repre- 
sentative of the "slowly" variable state (Zl) and added the flar- 
ing component (Z2) to match the April 1 1 and 16 high states. In 
this way the stability of the flux at energies <1 keV is easily ex- 
plained by the fact that the emission from Z2 is negligible when 



compared with that of Zl. The results of these calculations are 
represented in the SEDs of Fig. 11. We met some difficulty in 
obtaining a solution capable of modelling the observed spectral 
curvature mostly in the highest state. In fact, to deplete the SR 
from Z2 below ~1 keV, we need an electron population with 
a very low emissivity in this band. This can be obtained only 
by introducing a low energy cut-off, with the consequence of a 
high curvature of the low energy portion of IC spectra at TeV 
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energies. The values of parameters for this 2-zone SSC model 
are also given in Table 1. 

These results can be used to obtain some information on 
the acceleration process. According to Eq.(lO) we can derive 
the energy gain as a function of r. 

s = 10«/2'' (23) 

which, given that r ^ 1.0, gives s ^ 3.2^. This implies that, 
unless q is very small, the energy gain is about a factor 3. 
Consequently, the number of steps to increase the energy up 
to y ^ 10^ is around 10. The acceleration must therefore be 
efficient and should probably occur for rather short times. 

It is interesting to compare our SEDs of the 1997 flare of 
Mkn 501 with those calculated by other authors who adopted 
different physical models. Generally, there is no explicit eval- 
uation of a spectral curvature parameter for both SR and IC 
components and consequently the electron energy distribu- 
tions are assumed to be power laws with exponential cut-off 
(Konopelko et al. 2003, Krawczynski 2002) or a broken power 
law (Katarzynski 2001). 

To reproduce the spectral curvature in IC spectra, these au- 
thors need to introduce a heavy contribution from EBL inter- 
action with TeV photons. We stress that in our analysis the SR 
curvature is intrinsic, and we do not need break or cut-off to re- 
produce it. More interestingly, the IC curvature is also intrinsic 
and can be reproduced both using and neglecting the TeV pho- 
ton absorption by EBL photons. In particular in the, in case of 
large flares where the peaks of SR and IC components increase 
their energy, models with break or cut-off do not predict intrin- 
sically curved IC spectra, and so the spectral shape observed at 
TeV and the IC peak frequency are mainly modulated by what 
EBL model is chosen. Note also that EBL attenuation could 
consistently modulate the position of the IC peak because its 
optical depth varies with the energy and redshift. Observations 
of SR and IC components in objects with different redshifts and 
the study of the variations of their peak frequencies would be 
very useful for more accurate modelling of the EBL spectrum. 

It is interesting to compare the values of the main model 
parameters. As a further point Konopelko et al. (2003) and 
Krawczynski et al. (2002) used 6 values of 50 and 45, respec- 
tively, while Katarzynski et al. (2001) used 14, very close to 
our choice. 

6. Summary and conclusion 

Over the past decade, broad-band observations of blazars have 
clearly shown that their non-thermal SEDs are characterised 
by mildly curved shapes that can be approximated by simple 
power laws only in frequency intervals that generally not wider 
than a single decade. Starting from the first paper by Landau et 
al. (1986), we verified that a log-parabolic law is instead able 
to represent blazars' SEDs over several frequency decades. An 
advantage of this law is that it has only three free parameters. 
Furthermore, we have shown that for statistical acceleration 
mechanisms in which the acceleration probability is a decreas- 
ing function of the particle energy, as in Eq. (7), the resulting 
spectrum follows a log-parabolic law. Spectral properties of the 



SR can be derived easily only by using the <5-function approxi- 
mation. 

To find the exact relations between the spectral parame- 
ters of the energy distribution of the relativistic electrons and 
of their SR, we used a precise numerical code specifically de- 
veloped for blazar applications. The main results can be sum- 
marised as follows: 

- The SR spectrum is log-paraboUc with a curvature parame- 
ter approximately equal to the electron curvature parameter 
multiplied by 0.2; 

- spectral evolution due to SR cooling maintains a log- 
parabolic distribution whose curvature parameter increases 
linearly with time; 

- the IC SED around the peak can be approximately de- 
scribed by a log-parabola having a curvature parameter that 
ranges from values lower than the SR spectrum up to that of 
the electron distribution, depending upon whether the pho- 
ton scattering is in the Thomson or Klein-Nishina regime. 
The difference between the SR and IC curvatures can be 
used to discriminate between the two regimes. 

We applied our code to the calculation of the SR and IC 
spectra of Mkn 501 when it was simultaneously observed by 
BeppoSAX (Paper II) and CAT (Djannati-Atai et al. 1999) dur- 
ing the large outburst of April 1997. We estimated the spec- 
tral parameters of the emitting electrons from X-ray data and 
then computed the IC spectrum in the TeV band, taking the ab- 
sorption from pair production against EBL into account. Our 
results showed that the curvature at high energy is mostly in- 
trinsic and larger than in the X-ray range, which shows that 
the IC scattering is taking place in the Klein-Nishina regime. 
This finding suggests that intergalactic absorption of TeV pho- 
tons may be small, as shown by Dwek and Krennich (2005): 
their LLL spectrum could be a satisfactory description of EBL, 
but we cannot exclude that it has an even lower intensity. This 
interpretation agrees with a previous work by Krawczynski et 
al. (2000), where the low intensity DEBRA model (Malkan & 
Stecker 1998) was favored over another model of higher in- 
tensity. It is also supported by a recent paper by Aharonian et 
al. (2005a), who used H.E.S.S. data to find that in the time- 
averaged observations of Mkn 421, the cutoff energy of the IC 
spectrum is at ~3.1 TeV, lower than that found in Mkn 501 
(about 6.2 TeV). Considering the nearly equal redshifts of the 
these two blazars, they conclude that the cutoff is not due to 
EBL attenuation but is intrinsic to the sources, in agreement 
with our conclusion about the IC curvature. We conclude (i) 
that the extragalactic space is probably more transparent to TeV 
photons than previously assumed and (ii) that a larger number 
of blazars at higher redshifts than previously thought could be 
detected in this range. 

The same conclusion has been reached very recently by 
Aharonian et al. (2005b), who report the detection of the blazar 
lES 1101-232 at z = 0.186 (by H.E.S.S. in the TeV band). An 
important consequence of the intrinsic curvature of TeV spec- 
tra is that it can be variable in time and possibly related to the 
X-ray curvature. As a final remark, we stress the relevance of 
the simultaneous broad-band observations that are necessary 
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to obtain a reliable estimate of curvature parameters for both 
components. 
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